Spectral densities of superconducting qubits with environmental resonances. 
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In this paper we derive the environmental spectral den- 
sity for flux [l| , phase Q and charge Q qubits when each 
of them is coupled to an environment with a resonance. 
From the spectral density we obtain the characteristic 
spontaneous emission (relaxation) lifetimes T± for each 
of these qubits, and show that there is a substantial en- 
hancement of Ti beyond the resonant frequency of the en- 
vironment. The circuits considered are shown in Fig.[TJ[2] 
and [3] 

The flux qubit shown in Fig. [T] is measured by a dc- 
SQUID. Hence to study decoherence and relaxation time 
scales, one has to consider the noise that is transferred 
from the qubit to the dc-SQUID. 
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FIG. 1: Flux qubit measured by a dc-SQUID (gray line). 
The qubit is in the inner SQUID loop with critical current I c 
and capacitance C.j for both junctions. The inner SQUID is 
shunted by a capacitance C„, and environmental resistance R 
and is biased by a ramping current lb- The dc-SQUID loop 
has junction capatitance Co and critical current I c o- 

The classical equation of motion for the dc-SQUID is 
C 4>+ ^J c0 sin<£- 5^1+ / dt'Y(t-t')<j>(t')=0 (1) 
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where <j> is the gauge invariant phase across the Josephson 
junction of the outer dc-SQUID loop, I c q is the critical 
current of its junction, <£>o = h/2e is the flux quantum, 
and the total induced current in the outer dc-SQUID is 
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Here 4> p and <f> m are the sum and difference of the gauge 
invariant phases across the junctions of the inner SQUID, 
Ldc is the self- inductance of the inner SQUID, and J\ is 
the bilinear coupling between <fi m and <p p at the potential 
energy minimum. The term <5</>o = irMqI C ir/<&o, where 
I C ir is the circulating current of the localized states of the 



qubit (described in terms of Pauli matrix a z ), and M q is 
the mutual inductance between the qubit and the outer 
dc-SQUID. The last term in Eq. ([T]) is the dissipation 
term due to effective admittance Y(cu) felt by the outer 
dc-SQUID. 

For the charge qubit shown in Fig. [2] the classical equa- 
tion of motion for the charge Q is 

\ 1-uj^Lj(uj)Cj C g J 

Here, V g is the gate voltage, C g is the gate capacitance, 
Lj and Cj are the Josephson inductance and capacitance 
respectively, Z(u) is the effective impedance seen by the 
charge qubit due to a transmission line resonator (cavity) , 
f2 is the resonant frequency of the resonator, and Q is the 
charge across C g . 
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FIG. 2: Circuit diagram of the Cooper-pair box. The super- 
conducting island is connected to a large reservoir through 
a Josephson junction with Josephson energy Ej and capaci- 
tance C g . The voltage bias V g is provided through a resonator 
(cavity) having environmental impedance Z(lj), which is con- 



For the phase qubit shown in Fig. [3] the classical equa- 
tion of motion is 

Col+—I c0 sm 1 -—I+ / dt'Y(t-t')i(t') = (4) 

where I c q is the critical current of Josephson junction J 
in FigfTJ I is the bias current, and <&o = h/2e is the flux 
quantum. The last term is the dissipation term due to 
Y(u>) which is the effective admittance as seen by the 
dc-SQUID. 

The equations of motion described in Eqs. (HJ) , (J3J) ? and 
0, can be all approximately described by the effective 
spin-boson Hamiltonian 

H = ^cr z +J2^kblb k +H 3B , (5) 
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The spectral density for the charge qubit shown in 
Fig. [2] is obtained by solving for the normal modes of 
the resonator and transmission lines, including an input 
impedance R at each end of the resonator. It is given by 
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FIG. 3: Schematic drawing of the phase qubit with an RLC 
isolation circuit. 



written in terms of Pauli matrices Ui (with i = x,y, z) 
and boson operators bk and b\ . The first term in Eq. (J5]) 
represents a two-level approximation for the qubit (sys- 
tem) described by states |0) and |1) with energy differ- 
ence fiujQi . The second term corresponds to the isolation 
network (bath) represented by a bath of bosons, where 
bk and bl are the annihilation and creation operator of 
the fc-th bath mode with frequency u>k- The third term 
is the system-bath (SB) Hamiltonian which corresponds 
to the coupling between the environment and the qubit. 

At the charge degeneracy point for the charge qubit 
(gate charge N g = 0), at the flux degeneracy point (ex- 
ternal flux <f> oxt = 7T$ ) f° r the flux qubit, and the suit- 
able flux bias condition for the phase qubit (external flux 
$ a = Li^o) H S b reduces to, 



H SB = ^a x n{l\v\Q) 
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where v = <f> for the flux qubit, Q for the charge qubit, 
and 7 for the phase qubit. The spectral density of the 
bath modes J(w) = h J2k ( w — nas dimensions of 
energy and can be written as J(u>) = LuReY(uj)(&o/2TT) 2 
for flux and phase qubits and J(u>) = 2hcue 2 /KReZ(cu) 
for charge qubits. 

For the flux qubit circuit shown in Fig. [TJ the shunt 
capacitance C s is used to control the environment, while 
the Ohmic resistance of the circuit is modelled by R. In 
this case, the environmental spectral density is 
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when u> m ^> max(w p , loq), and when the dc-SQUID is far 
away from the switching point to be modelled by an ideal 
inductance Lj. Here, fii = \/\fLjC s = y/iirl c / (C s $o) 
is the plasma frequency of the inner SQUID and plays 
the role of the resonant frequency, where I c is the crit- 
ical current for each of two Josephson junctions. Also, 
Ti = l/(C s R) plays the role of the resonance width, and 
a x = 2(eI p I b M q ) 2 /(C 2 h 2 Rnf) reflects the low frequency 
behavior. The coupling between the flux qubit and the 
outer dc-SQUID occurs emerges from the interaction of 
the persistent current I p of the qubit and the bias current 
h of the dc-SQUID via their mutual inductance M q . 



were is the resonator frequency, I is resonator length, 
c is the capacitance per unit length of the transmission 
line, C g is the gate capacitance, and Cj is the junction 
capacitance. The quantity T 2 = Sl 2 /Q where Q is the 
quality factor of the cavity. 

For negligeable Y int (uj), the spectral density of the 
phase qubit shown in Fig. (|3|) is 
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where a 3 = L 2 /((L + Li) 2 R) w (L/L^/R is the 
leading order term in the low frequency ohmic regime, 
3 = y/(L + Li)/(LLiC) w 1/VLC is essentially the 
resonance frequency, and T 3 = \/{2CR) plays the role of 
resonance width. Here, we used L\ L corresponding 
to the relevant experimental regime. 

Once the spectral functions of the environments are 
known, the relaxation rates are derived following stan- 
dard methods At finite temperatures, the relaxation 
time 



Ti, 



= j3iJi(ijj) coth 
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is directly related to the environmental spectral den- 
sity. Here, the index i — 1,2,3 labels the results for 
flux, charge and phase qubits, respectively The param- 
eter f3i takes values /3\ = 1 for the flux qubit of Fig. [TJ 
P2 = C 2 /{C g + Cj) 2 for the charge qubit of Fig. [2] and 

/3 3 = 1/ ($ /27t) 2 C w i for the phase qubit of Fig. [3] 

An inspection of Eq. (|10p and the corresponding spectral 
functions Ji{uj) shows that suitably detuning the qubit to 
higher frequencies beyond the environmental resonance 
can enhance the characteristic spontaneous emission life- 
time (relaxation time) of the qubit by a couple of 
orders of magnitude as discussed in Ref. [5j. 
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